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INTRODUCTION 


XYZPLOT  is  a  passive  graphics  program  for  displaying  the  results  of 
potential  flow  calculations.  It  was  written  for  the  CDC  6000  series 
computers  and  the  Stromberg-Carlson  4060  Stored  Program  Recording  System. 
XYZPLOT  was  specifically  programned  to  mate  with  the  input  formats  and 
restart  files  of  the  XYZ  Potential  Flow  Program  (XYZ  PFP),1  but  with  proper 
interfacing  XYZPLOT  may  be  used  to  display  results  from  similar  programs. 

The  body  may  be  specified  in  surface  element  form  as  in  XYZ  PFP,  or 
in  an  alternate  form  (outline  form)  in  which  the  user  specifies  points 
which  comprise  an  outline  of  the  body.  However,  when  plotting  isobars  or 
on-body  streamlines,  the  surface  element  points  must  always  be  specified. 

A  Mach  nunber  correction  is  included  for  plotting  isobars  in  the 
case  of  subsonic  flow  in  air.  A  special  version  of  the  XYZ  PFP  has  been 
developed  to  generate  pressure  coefficients  for  this  case.* 

XYZPLOT  is  intended  to  supersede  two  earlier  programs,  VIEKPLOT  and 
2 

ISOPLOT.  It  has  been  successfully  applied  to  several  problems,  including 
a  study  of  isobars  on  helicopter  rotor  hubs  and  a  comparison  of  the 
XYZ  PFP  with  analytic  results  for  an  ellipsoid. 

XYZPLOT  and  its  interactive  graphics  counterpart,  DESIGN, 3,+  represent 
the  present  graphics  capability  for  the  XYZ  Potential  Flow  Program. 


Consult  Code  1843,  Numerical  Fluid  Dynamics  Branch,  for  details. 

+  DESIGN  presently  has  the  capability  to  display  and  edit  body  points. 
Streamline  and  isobar  routines  will  soon  be  added. 

*  Dawson,  Charles  W.  and  Dean,  Janet  S. ,  "The  XYZ  Potential  Flow  Program," 
Naval  Ship  Research  and  Development  Center  Report  3892  (June  1972) . 

2 

Reingruber,  John  K. ,  "Three-Dimensional  Computer  Graphics  for  Fluid 
Dynamics,"  Confutation  and  Mathematics  Department  Technical  Note  CMD-56-72 
(November  1972) . 

3  Kelly,  Barbara  M.  and  Marquardt,  Mary  Beth,  "Interactive  Helicopter 
Design:  Geometry  Package  User's  Manual,"  Computation  and  Mathematics 
Department,  Departmental  Report  CMD-28-74  (September  1974). 
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DESCRIPTION  OF  METHODS 


PROJECTION 

XYZPLCTT  projects  points  in  space  onto  a  plotting  plane,  S,  by  one 
of  two  user  selected  methods.  Both  methods  require  the  specification  of 
on  observation  point,  Qq  ■  (x0,y0,zQ),  as  well  as  two  points, 

Ql  *  and  Q2  -  (x2,y2,z2),  which  lie  in  the  plotting  plane.  Qj 

is  taken  as  the  origin  of  the  plotting  plane  coordinate  system 
(designated  by  xs,ys),  and  Q2  determines  the  positive  xg  direction.  A  point 
lying  in  the  positive  y  direction  Q,  -  (x,,y,,z.)  is  calculated  by  finding 
the  vector  product  of  Q.Qq  and  QjQ,.  No  distortion  will  occur  when 
is  perpendicular  to 

Perspective  Projection 

The  perspective  projection  Qp  -  (x^,yp,Zp)  of  a  point  in  space 
Qg  ■  (Xg./g.Zg)  is  the  intersection  of  the  line -of -sight  (QqQ^)  wi th  the 
plotting  plane.  See  Figure  1.  This  may  be  formulated  as  follows: 


^  i 


Cl) 

(2) 


Constraint  (1)  is  satisfied  if 


(--<  t  <«) 


In  Cartesian  coordinates, 


y*p  ■ 

Vyp  ■  t(yo-V 

Vzp  ■  t(W 


■  m 


Figure  1.  Geometry  for  Perspective  Projection 

The  plotting  plane,  S,  contains  points  Qj.Q.,,  anJ  which 
define  xs,yg  coordinate  system.  The  observation  point  Qq 
and  a  point  in  space  Qg  determine  the  line-of-sight .  The 
projected  point  Qp  lies  at  the  intersection  of  the  plotting 
plane  with  the  line-of-sight. 


or 

t .  3Li  -  r*ZE .  !£!e 
Vxg  yo'yg  V:g 

This  relation  yields 

(y0-yg>txg'V  ■  WY'V 
(xg'V(Io-V  *  (vxg><yV 

(YyP)(VV  •  <Vyg>(W 

which  may  be  rewritten  as 


lyo-yg)xP  * 

(VVyP 

"  y0xg  -  Vg 

(3) 

*  (Yxo,2P 

*  z0xg  ‘  Vg 

(4) 

(Y  2o)yP  *  (yo-yg)2p 

•  yo2g  '  vg 

(5) 

Only  two  of  these  equations  are  needed  to  define  the  line;  the  third 
equation  is  always  derivable  from  the  other  two. 

Constraint  (2)  is  satisfied  if 

■  o 

In  Cartesian  coordinates, 

(yxiUxo-Xj)  -  (yp-yjHyQ-yi)  -  (yV(vV  *  0 
or 

(Vxi)xp  ♦  (Xo-y^Xp  -  (yyzp  -  (x0-x1)x1 

16) 

4  (y0-yi)yi  4  (v2i)zi 

Equation  (b)  along  with  two  equations  selected  from  equations  (3), 
(4),  and  (5)  constitute  a  system  of  three  equations  in  the  three 
unknowns  Xp,  yp,  and  xp.  The  system  will  have  a  solution  provided  ^ 
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is  distinct  from  Q  and  is  not  parallel  to  the  plane  S. 


'R  'U'g 

Orthographic  Projection 

The  orthographic  projection  =  (x  .y^.y  of  a  point  in  space 


Q  *  (x  ,y  j  is  the  projection  of  Q  onto  the  plotting  plane  along  a 
K  R  8  K  8 

line  normal  to  the  plane.  See  Figure  2.  Tins  may  he  formulated  as 


*5?  ■  5??  • 11 

■  5??-  0 

5C?  '  5??  ’  0 

(for  orthographic  projection  is  taken  as  normal  to  SI. 
In  Cartes ian  coordinates, 

lyxpHvV  +  (>'g'>pu>’:'>'i)  +  [\-zrnzyzi]  = 

(yxp)(x3-xi)  ♦  (>yyp H/j-yp  ♦  = 


VxiUW  -  l>p-yi)(y0->'i)  -  = 


which  may  be  rewritten  as 


(x2'xi'V (1'2'i'1*V<!J’:i1:p  ’  lx:'xi)V(>Vyi)Vu: 


:i'=e 


(Xj-xpx  H/j-ypy  'Uj-zpr  ■  (X.-X,)X  -O-j-ygy  .(--3 


:i':y 


(x0-xl)V(y0->,l)VU0-‘l)2p  ■  lx0'xl)V(Vyl,yl*u0' 


=l’:l 


or 


ix2-xi^  (y 2 i ^  (z2~zi 

(Xj-Xj)  (yj-yj)  (:3'zi 

cIXq-Xj)  (y0->’il  C-0" zi 


ix^xpxgHyyypypn 

OvVyOy^Y1 

ixo-xpXiHy0-yi)yin 


-  i  - 
“i  g' 


V:1):C 


"0  “1  1> 


This  system  of  linear  equations  has  no  solution  if  the  coefficient  matrix 
is  singular,  i.e.,  when  •  (Q-0^  X  QpQ^l  =  0.  There  is  no  d.uigi 

of  this  happening  as  long  as  is  nearly  normal  to  S. 
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Figure  2.  Geometry  for  Orthographic  Projection 

The  plotting  plane,  S,  contains  points  Qj,Q7,  and  which 
define  the  xs,ys  coordinate  system.  The  observation  point 
is  no  al  to  S  at  A  point  in  space  Q  is  projected  along 

a  line  normal  to  S.  The  projected  point  Q  ^  lies  at  the 
intersection  of  this  normal  with  S. 


Conversion  to  x  ,y  Coordinates 
s  '  s 


Cnee  the  three  coordinates  of  the  projected  point  Qp  are  determined 
(by  either  projection  method)  it  is  necessary  to  find  the  corresponding 
coordinates  in  the  x$,ys  coordinate  system.  This  is  done  by  taking  the 
scalar  product  of  the  position  vector  in  the  x 


respective  base  vectors  xg  -  and  yg 

The  coordinates  are  then  normalized  by  dividing  by  R,  where 


C* 


yg  system  (Q^)  with  the 


Q3  is  chosen  so  that 


i*,i  “  ly.i 


IxJ  -  lye 


*p. 


-  K-^ - 


— — 


See  Figure  3. 


F’gure  3.  Conversion  to  x$ ,y$  Coordinates 


LOCATING  HIDDEN  SURFACE  ELB4EMTS 


Associated  with  each  surface  element  is  a  normal  Ft  specified  at  the 
geometric  center  C  of  the  element.  Ft  is  defined  by  the  vector  product  of 
the  diagonals  of  the  surface  element.  See  Figure  4. 


i+lj 


i+lfj+l 


PiJ+l 


Figure  4.  General  Surface  Element 


N  =  (Pi,j  Pi+l,j+l)  x  (Pi,j+l  Vl.j3 

The  coordinates  of  C  are  found  by  taking  the  average  of  the  coordinates 
of  the  surface  element's  comer  points. 

YC  ‘  (yi , j , j ♦l*yi*l , j *yi*l 

Neither  of  these  calculations  is  required  when  XYZPLOT  is  used  to 
plot  isobars  because  the  normals  and  centroids  are  available  on  a  disk  file 
generated  by  the  XYZ  PFP. 

If  the  normal  of  a  given  quadrilateral  points  away  from  the  observer, 
it  is  considered  hidden.  This  is  determined  by  finding  the  scalar  product 
of  the  normal  with  the  vector  indicating  the  direction  of  the  observer,  D. 

SP  ■  ft  •  if 
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During  perspective  projection  D  is  the  vector  from  the  observation 
point  to  the  center  of  the  surface  element.  During  orthographic  projection 
it  is  the  normal  to  the  plotting  plane.  If  SP>0,  the  element  is  hidden  and 
therefore  not  plotted.  This  method  does  not  detect  surfaces  which  are 
hidden  due  to  blockage. 

DETERMINING  HIDDEN  STREAMLINES 

Hidden  off -body  streamlines  are  determined  as  follows:  The  user 
specifies  a  plane  by  its  origin  and  a  normal  which  points  toward  the 
observer.  Any  streamline  segment  having  an  endpoint  which  is  not  on  the 
observer  side  of  the  plane  is  not  plotted.  This  leaves  the  user  the  option 
of  selecting  which  streamlines  will  be  plotted  in  each  view. 

XYZPLOT  automatically  suppresses  on-body  streamlines  which  lie  on 
hidden  surface  elements. 


PLOTTING  ISOBARS 

Determining  Pressure  Coefficients 

The  local  velocities  (at  the  null  points)  are  calculated  from  a  given 

onset  velocity  ( VB  ,  ,  VB  )  and  the  local  velocities  for  the  X,  Y,  and 

°°x  "y  “z 

Z  flows  ((VXl,VYl,VZl)i,  (VX2,VY2,VZ2)i,  (VM  ,VY3  ,VZ3)  .  The  local 

velocities  for  the  X,  Y,  and  Z  flows  are  written  on  TAPE03  by  section  4  of 
the  XYZ  PFP.  At  each  null  point,  i,  the  local  velocities  are 


From  these  local  velocities,  the  pressure  coefficients  at  each  null  point 
a.  e then  calculated.  For  inconpressible  flow 


1 


VX?  +  VY?  +  VZ? 

?  +  v2  TV 


x  y  z 


(7) 
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If  the  compressible  flow  version  of  the  XYZ  PFP  was  used,  a  Mach 
nunber  correction  is  applied 


Locating  Isobars 

Associated  with  each  null  point  is  a  pressure  coefficient  (computed 
from  either  liquation  7  or  liquation  8).  The  set  of  null  points  (and  there¬ 
fore  the  set  of  pressure  coefficients)  form  a  system  of  quadrilaterals.  This 
system  is  shown  as  broken  lines  on  the  body  section  in  Figure  5. 


Figure  5.  Pressure  Coefficient  Quadrilateral  System 
A  specified  isobar  (Cp  )  is  located  by  searching  over  the  vertices 

*3 

of  each  of  these  quadrilaterals  to  find  cases  where  the  comer  values 
bracket  Cp  .  Once  found,  linear  interpolation  will  yield  at  least  two 
points  on  that  quadrilateral  having  approximately  the  specified  pressure 
coefficient.  These  points  are  connected  to  form  an  isobar  segment.  If 
more  than  two  points  along  the  quadrilateral's  edge  are  fi  md  to  have  the 
same  pressure  coefficient,  the  pairs  having  the  smallest  combined  lengths 
are  arbitrarily  connected.  See  Figure  6. 
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Figure  6.  Hypothetical  Pressure  Coefficient  Grid  with  Isobars 

The  scheme  used  to  index  the  surface  elements  (and  null  points)  does 
not  carry  across  section  boundaries,  therefore  preventing  isobars  from  being 
drawn  between  two  sections.  A  simple  option  is  available  to  the  user  which 
often  circumvents  this  problem.  Consider  the  situation  represented  in 
Figure  7. 


Figure  7.  Sectional  Interface 


It  appears  that  the  isobar  ending  at  point  C  should  continue  across 
the  section  boundary  to  point  D.  With  the  connection  option  on,  point  B 
will  be  connected  to  point  E  provided  that 

|ri|  ♦  \&\  >  \tf\ 

If  more  than  one  point  in  section  2  meets  this  requirement,  the  nearest  point 
(to  C)  is  used. 

This  method  will  never  connect  two  points  representing  isobars  of 
different  value,  but  it  may  miss  some  appropriate  connections  or  improperly 
connect  different  isobars  of  the  same  value. 

PREPARATION  OF  INPUT 

BODY  DESCRIPTION 

XYZPLOT  will  plot  a  body  specified  by  either  surface  elements  or 
line  segments.  However  the  following  restrictions  arise  whenever  on-body 
streamlines  or  isobars  are  to  be  plotted:  1)  the  body  must  be  plotted  in 
surface  element  form  when  on-body  streamlines  are  being  plotted;  21  the 
surface  element  points  must  be  read  in  when  isobars  are  going  to  be  plotted, 
even  if  outlining  is  used  for  body  representation.  The  user  may  insure 
compliance  with  these  restrictions  by  checking  his  final  selection  of 
parameters  with  the  valid  combinations  shown  in  Table  1. 

Surface  Elements 

The  surface  elements  are  described  in  the  same  manner  as  in  the  XY2  PFP, 

4 

as  explained  by  the  following  excerpt  from  Dawson  and  Dean. 

"The  surface  of  the  body  is  divided  into  sections.  Each  section  is 
described  separately  from  the  other  sections  so  that  the  description  of  one 
section  can  be  changed  without  requiring  changes  in  the  descriptions  of  the 
other  sections.  Standard  sections  can  thus  be  described  and  various  combina¬ 
tions  cf  these  sections  can  be  used  in  many  problems  with  changes  introduced 
as  needed. 

4 

Dawson  and  Dean,  op.cit.,  pp.  30-33. 
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TABLE  1  -  VALID  PARAMETER  COMBINATIONS 


Combination 

ITYPE 

ISE 

IOL 

I  BODY 

ISTRM 

1 

0 

0 

4,5 

2 

0 

2 

0 

2,5 

0 

1 

0 

3 

0 

2,5 

4,5 

1,2 

0 

4 

1 

0 

4,5 

2 

0,1,2 

5 

1 

2,5 

0 

1 

0,1,2 

6 

1 

2,5 

4,5 

1,2 

0,1,2 

7 

2 

0 

4,5 

2 

0 

8 

2 

2,5 

0 

1 

0,1,2 

9 

2 

2,5 

4,5 

1 

0,1,2 

10 

2 

2,5 

4,5 

2 

0 

11 

3 

0 

4,5 

2 

0 

12 

3 

2,5 

0 

1 

0,1,2 

13 

3 

2,5 

4,5 

1,2 

0,1,2 

"Each  section  is  subdivided  by  two  sets  of  lines,  more  or  less  at 
right  angles  to  each  other.  The  input  consists  of  the  X,  Y,  and  Z 
coordinates  of  the  points  at  the  intersection  of  the  lines  and  the  indices 
N,  M  to  indicate  the  lines.  The  indices  for  the  lines  should  be  chosen 
so  that  when  viewed  from  outside  the  body,  N  increases  to  the  right  when  V 
increases  toward  the  top.  Any  rotation  of  this  arrangement  is  also 
permitted  but  a  reflection  will  result  in  an  inward  pointing  normal  vectoi 
Of  course,  if  a  left-handed  coordinate  system  is  being  used,  N  and  M 
must  be  interchanged...  .  [See  Figures  8,  9,  and  10.] 

"The  plane  quadrilaterals  which  make  up  the  approximation  to  the 
surface  are  generated  from  the  point  data.  A  quadrilateral  will  be 
generated  only  when  the  four  comer  points  N,M;  N,M+1;  N+l.M+l;  and  N*1,M 
are  given.  Thus,  a  section  with  a  hole  in  it  or  a  section  with  two  parts 
can  be  specified.  The  advantage  of  this  arrangement  is  that  part  of  a 
body  can  be  changed  or  more  accurately  represented  without  requiring  that 
all  the  points  be  renumbered." 
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SECTION  2 


Figure  10  Representation  of  a  Submarine 


(Xitline 

The  body  outline  is  divided  into  components,  and  each  component  is 
composed  of  line  segments.  The  line  segments  are  specified  by  two 
consecutively  indexed  points.  To  plot  the  outline,  a  line  is  drawn  from 
the  point  indexed  M  to  the  point  indexed  M+l  within  a  component.  To  make 
a  component  close  on  itself,  the  point  with  the  largest  index  must  be 
identical  to  the  point  with  the  smallest  index. 

DATA  ORGANIZATION  AND  FORMAT 

XYZPLCTT  uses  a  combination  of  cards,  disk  files,  and  tapes  as  input. 

It  requires  a  problem  identification  card,  a  control  parameter  card,  body 
point  cards,  and  a  user  identification  card.  The  body  point  cards  mav  be 
punched  cards,  or  card  images  which  reside  on  a  disk  or  tape  file,  Each 
view  then  requires  a  set  of  plotting  information  cards  followed  bv  an 
end-of-record  card.  Streamline  or  isobar  data  are  read  from  a  disk  or  tape 
file. 

Problem  Identification  and  Control  Parameters 

The  problem  identification  card  is  first  in  the  data  deck  and  is  m 
FORTRAN  format  (6A10) .  The  information  in  this  card  should  be  approx i mat el v 
centered  between  columns  2  and  60.  The  control  parameter  card  is  next  and 
is  in  format  (815,  F5.3).  All  parameters  must  be  right -justified,  as  blanks 
are  interpreted  as  zeros.  This  card  must  contain  the  following  information: 


PARAMETER 

COLUMN 

DESCRIPTION 

NVIEWS 

1-5 

Number  of  views  to  be  gene  rat  oil. 

ISE 

10 

Surface  element  tape  indicator 

ISE 

SITUATION 

0 

If  not  reading  surface  elements. 

y 

4 

If  reading  from  disk  file  or 

tape  (attach  as  TAPEZj . 

5 

If  reading  from  cards. 
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PARAMETER 

IOL 

ITYPF 


I  SYS 

IVERS 
I  OUT 

NQE 

XMAC 


COLIMJ  DESCRIPTION 

15  Oat line  tape  indicator. 

IOL  SITUATION 

0  If  not  reading  outline  points. 

4  If  reading  from  disk  file  or  tape 
(attach  as  TAPE4) . 

5  If  reading  from  cards. 

20  Flow  data  type  indicator. 

ITYPE  SITUATION 

0  No  data  (only  body  may  be  plotted) . 

1  Off-body  streamline  points  (e.g. 
TAPE16  from  section  6  of  XYZ  PFP) . 

2  On-body  streamline  points  (e.g. 
TAPE17  from  section  7  of  XYZ  PFP) . 

3  Isobar  Information  (e.g.  TAPE03 
from  section  4  of  XYZ  PFP) . 

25  Coordinate  System  orientation 

indicator. 

I SYS  SITUATION 

-1  Left-handed. 

1  Right-handed. 

30  Version  of  XYZ  PFP  used. 

35  Set  IOUT=l  to  prevent  input  echo- 

checking. 

36-40  Number  of  surface  elements.  (Set 

NQE=0  if  not  reading  surface 
element  points). 


41-45  Mach  number.  Set  XMAOO  unless 

compressible  flow  version  of 
XYZ  PFP  was  used. 
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Body  Point  Cards 

It  is  possible  to  have  the  body  plotted  in  surface  element  form  in 
some  views  and  in  outline  form  in  other  views  during  a  given  computer  run. 

To  do  this,  surface  element  points  and  outline  points  must  both  be  read  in 
(in  that  order) .  The  user  must  be  careful  to  specify  the  proper  values 
for  ISE  and  IOL  on  the  control  parameter  card. 

Surface  Element  Points.  If  the  body  is  to  be  plotted  in  surface 
element  form  in  any  view,  or  if  isobars  or  on-body  streamlines  will  be 
plotted  in  any  view,  the  surface  element  point  cards  must  be  next  in  the  data 
deck.  These  cards  are  the  same  as  the  body  point  cards  for  the  XYZ  PFP, 
except  that  the  last  parameter  on  each  card  (the  normal  component  of  the 
velocity  at  the  body  surface)  is  not  required.  The  format  for  these  cards 
is  (3F12.9,  414).  The  maximum  number  of  body  sections  allowed  is  25,  the 
maximum  number  of  body  points  allowed  in  a  body  section  is  726,  and  the 
entire  body  may  consist  of  up  to  650  surface  elements.  Each  point  card  must 
contain  the  following  information: 


PARAMETER 

COLUMN 

DESCRIPTION 

XI 

1-12 

X-coordinate  of  the  point. 

YI 

13-24 

Y-coordinate  of  the  point. 

ZI 

25-36 

Z-coordinate  of  the  point. 

NI 

39-40 

N-index  of  the  point  (1  ?  NI  5  70) . 

MI 

43-44 

M- index  of  the  point  (1  -  MI  5  40). 

NS 

45-48 

Section  identification  number  (1‘ NS' 9999) 
NS  must  be  unique  for  each  section 

NE 

52 

Change  indicator  for  the  direction  of  the 
normal  vector.  When  NE  on  the  first  caid 
of  a  section  is  not  blank  or  zero,  NI  and 
MI  are  interchanged  for  that  section  to 
change  the  direction  of  the  normal  vecto* 
NE  is  ignored  on  the  other  cards. 

All  point  cards  in  a  section  must  be  placed  together,  but  they  may  be  in  any 
order  within  the  section.  A  blank  card  or  an  end-of-recoid  must  fellow  the 
last  section. 

If  the  surface  element  points  are  to  be  read  in  from  a  disk  file  or  tape, 
the  data  file  must  be  attached  as  TAPE2. 
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Outline  Point  Cards.  If  the  body  is  to  be  plotted  in  outline  form  in 
any  view,  the  outline  point  cards  must  be  next  in  the  data  deck.  The 
format  for  the  outline  cards  is  (3F12.9,  214).  The  maximum  number  of  outline 
points  allowed  in  a  body  component  is  100,  and  the  body  may  consist  of  up 
to  25  components .  Each  point  card  must  contain  the  following  information: 


PARAMETER 

COLUW 

DESCRIPTION 

XI 

1-12 

X-coordinate  of  the  point. 

YI 

13-24 

Y -coordinate  of  the  point. 

ZI 

25-36 

Z -coordinate  of  the  point. 

MI 

39-40 

Index  of  the  point  (MI  <  100) . 

NC 

41-44 

Component  identification  number 
(1  5  NC  f  9999) .  NC  must  be 
unique  for  each  component . 

All  point  cards  for  a  component  must  be  placed  together,  but  they  may 
be  in  any  order  within  the  component.  A  blank  card  or  an  end-of-record  must 
be  placed  after  the  last  component. 

If  the  body  outline  points  are  to  be  read  from  a  disk  file  or  tape, 
the  data  file  must  be  attached  as  TAPE4. 

Streamline  Points  and  Isobar  Data 

During  a  given  computer  run  XYZPLOT  can  plot  either  on -body  streamlines, 
off-body  streamlines,  or  isobars.  The  user  must  be  careful  that  the  type 
of  data  file  attached  (as  TAPE50)  agrees  with  the  value  of  ITYPE  on  the 
control  parameter  card.  The  files  must  be  written  unformatted  as  follows 
(all  variable  types  are  standard  FORTRAN  default) : 

Off -Body  Streamline  Points.  Each  set  of  off -body  streamlines  requires  * 
one  record  with  the  following  parameters: 

DESCRIPTION 

Number  of  off -body  streamlines  in  the 
set  (NOBP  5  726) . 

Number  of  time  steps  in  this  set. 

Indicator  that  signals  the  last  set 
of  streamlines:  IEND=0  when  another  set 
follows;  IENI>1  to  signal  the  last  set. 

X,  Y,  and  Z  components  of  the  onset 
velocity  for  this  set 


PARAMETER 

NOBP 

NST 
I  END 


VX 

VY 

vz 


The  streamline  points  are  -ontained  in  the  next  NST+1  records.  This 
pattern  continues  until  all  the  streamline  sets  are  on  the  file.  As 
indicated  above,  the  last  set  is  flagged  by  IENI>1. 

The  following  example  shows  how  the  file  would  be  written  for  the 
sample  set  of  streamlines  in  Figure  11.  Streamlines  formed  by  negative 
time  steps  are  treated  as  a  separate  set  from  those  formed  b  positive  time 
steps,  even  though  they  may  have  the  same  step  0  starting  points. 

Record  #1  Record  #2 


NOBP, NST,1 END, VX,VY,VZ 
For  First  Set 

(For  Positive  Time  Steps) [X ,Y ,Z] A  Q , ••• [X,Y,Z]£  Q, 

Record  #3  Record  #4 

[x,y,z]a_1,---,[x,y,z]1._1 

[X,Y,Z] A  2 , • • • , [X,Y,Z]E  £ 

Record  #5  Record  #6 

NOBP ,NST , IEND.VX ,VY ,VZ 
For  Second  Set 

(For  Negative  Time  Steps) [X,Y,ZJA  Q, •••, (X,Y,Z)l  0, 

- 

Record  #7  Record  #8 

[x,y,z]a.1,-*-,[x,y,z]e.1 

[X,Y,Z)A  2>  "  *  > [X,Y,Z] g_2 

END'  OF  FILE - [ 


This  file  may  contain  streamline  sets  for  several  onset  flows:  m  a 
given  view  XYZPLOT  will  search  for  and  plot  all  off-body  streamlines  available 
for  a  specified  onset  flow. 

This  file  is  automatically  written  as  TAPE16  by  section  6  of  XY2  PFP. 

On- Body  Streamline  Points.  Each  set  of  on-body  streamlines  requires 
one  record  with  the  following  parameters: 

PARAMETER  DESCRIPTION 

NLIN  Number  of  streamlines  in  the  set. 

vx 

^  X,  Y,  and  Z  components  of  the  onset 

VZ  velocity  for  the  set. 

This  record  is  followed  by  NLIN  records  which  contain  the  streamline  points 


These  records  are  organized  as  follows: 

NSLP, (X(I) ,Y(I) ,Z(I) ,NQUAD(I) ,I-1,NSLP) 


where  NSLP  is  the  number  of  streamline  points;  X,  Y,  and  Z  are  their 
coordinates;  and  NQUAD  is  the  number  of  the  surface  element  on  which  the 
point  lies. 
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STEP  NUMBERS 


Figure  11.  Sample  Set  of  Streamlines 


The  following  example  shows  how  the  file  would  be  written  for  the  sample 
set  of  streamlines  in  Figure  11.  For  this  example  streamlines  formed  by 
negative  time  steps  are  treated  as  a  separate  set  from  those  formed  by  positive 
time  steps. 


Record  #1 

Record  #2 

NLIN,VX,VY,VZ 
for  first  set 

NSLP ,[ X , Y , Z , NQUAD ]*. Q ,[ X , Y , Z , NQUAD ] ^ x ,[ X , Y , Z , NQUAD ] ^_ 2 

Record  #3 

NSLP,  [X,Y,Z,NQUAD]g_0,  [X.Y.Z.NQUADlg^,  [X,Y,Z,NQUAD]g_2 

Records  #4-6 

Record  #7 

_ rv  v  7  NOUAm  + 

NLIN,VX,VY,VZ 
for  second  set 

[A,  I  9L  ,1NI^U/VDJ  g_2 
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This  file  may  contain  streamline  sets  for  several  onset  flows:  in  a 
given  view  XYZPLOT  will  search  for  and  plot  all  on-body  streamlines 
available  for  a  specified  onset  flow. 

This  file  is  automatically  written  as  TAPE17  by  section  7  of  XYZ  PFP. 

Isobar  Data.  To  locate  and  plot  isobars,  XYZPLOT  requires  the  local 
velocities  at  the  null  points  for  the  X,  Y,  and  Z  flows.  These  are  read 
from  TAPE03  of  section  4  of  the  XYZ  PFP. 


User  Identification  Card 

The  next  card  is  the  user  identification  card.  It  is  in  format  ( 3A1 0 ) 


and  must  have  the 

following  information: 

PARAMETER 

COLUNM 

DESCRIPTION 

NAME 

1-10 

User's  last  name. 

CODE 

11-20 

User's  code. 

PHONE 

21-30 

User's  phone  number. 

Plotting  Information  Cards 

The  plotting  information  cards  are 

next,  and  must  be  provided  foi  cath 

view.  Each  set  must  be  followed  by  a  7/8/9  (end -of -record)  card.  The  f,vst 
six  plotting  information  cards  are  in  format  (6A10/4I5/3(3F20.7/) ,5F15.” ' . 

They  contain  the 

following  information  in  the  given  order: 

PARAMETER 

CARD/COLUMNS 

DESCRIPTION 

TITLE 

1/2-60 

Alphanumeric  field  containing  the 
title  to  be  displayed  on  the  frame 

I  BODY 

2/5 

Body  plotting  indicator. 

I BODY  SITUATION 

1  Body  will  be  plotted  in  surface 
element  form. 
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PARAMETER 


CARD/COLlftlS 


ISTRM 


I  GRID 


IMODE 


DESCRIPTION 

SITUATION 


I BODY  (cont.)  SITUATION 

2  Body  will  be  plotted  in  outline 
form. 


ISTRM 


IGRID 


IMODE 


3/1-20 

3/21-40 

3/41-60 

4/1-20 

4/21-40 

4/41-60 

5/1-20 

5/21-40 

5/41-60 

6/1-15 

6/16-30 

6/31-45 

6/46-60 


Streamline/ Isobar  indicator. 

SITUATION 

Plot  no  streamlines  or  isobars. 

Plot  streamlines  or  isobars. 

Advance  frame,  then  plot  stream¬ 
lines  or  isobars. 

Grid  Indicator. 

SITUATION 

No  grid. 

Grid  displayed. 

Advance  frame;  then  draw  grid. 

Mode  Indicator. 

MODE 

Perspective 

Orthographic 

Coordinates  of  the  observation  point. 

Coordinates  of  the  plotting  plane 
origin. 

Coordinates  of  a  point  along  the 
xg-axis  of  the  plotting  plane. 


Values  in  the  x_ 


coordinate 


system  of  the  left,  right,  bottom,  and 
top  margins,  respectively,  of  the 
plotting  frame.  If  not  specified, 
these  quantities  will  be  computed 
according  to  the  value  of  SR  below. 
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PARAMETER 

SR 


CARD/C0LIM4S 

6/61-70 


SR 

1.0 

0  <  SR  <  1.0 


0.0 


DESCRIPTION 

Size  ratio  which  governs  the  size 
of  the  picture. 

SIZE 

Full  size 

Fraction  of  full  size,  (e.g., 
if  SR-,5,  then  the  picture 
will  be  scaled  to  half-size). 
XL,  XR,  YB,  YT  are  specified 
above. 


IGRID  allows  the  user  to  select  a  rectangular  grid  to  appear  on  one  of 
the  output  frames.  The  grid  spacing  is  approximately  7  millimeters.  The 
parameters  ISTRM  and  IGRID  permit  each  view  to  occupy  up  to  three  frames  of 
output.  Table  2  may  be  used  to  predict  the  output  resulting  from  various 
combinations  of  ISTRM  and  IGRID. 
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TABLE  2  -  ISTRM  AND  IGRID  COMBINATIONS 


IF  ISTRM  - 

0 

0 

0 

1 

1 

1 

2 

9 

9 

AND  IGRID  - 

0 

1 

2 

0 

1 

2 

0 

1 

*> 

The  body  appears 
in  frame 

1 

1 

1 

1 

1 

1 

1 

l 

i 

The  streamlines  or 
isobars  appear  in 
frame 

1 

1 

1 

2 

2 

The  grid  appears 
in  frame 

1 

2 

1 

2 

i 

9 

The  seventh  plotting  informati 
the  following  information: 


PARAMETER 


COLUMJ 


card  is  in  format  (215)  and  requires 


DESCRIPTION 


NSP 

(NCP) 

4-5 

/  Number  of  body  sections  (components)  to 
/  be  plotted.  NSP  (NCP)  <  25. 

I»l 

10 

/  Number  of  planes  of  symnetry 

This  is  followed  wit 
PARAMETER 
NSN 


SYMMETRY  PLANES 

None 

Y  =  0 

Y  =  0  and  Z  =  0 

Y  =  0,  Z  =  0,  and  X  =  0 

the  following  card(s)  in  format  (1615) : 

DESCRIPTION 

Array  that  contains  the  section  (component) 
numbers  of  the  sections  (components)  to 
be  plotted.  The  dimension  of  this  array 
is  25. 

When  plotting  /sobars  or  streamlines,  the  following  cards  must  be 
included  when  appl i/able. 

The  ONSET  FU/)W  card  contains  the  onset  velocity  components  for  this 
view  in  format  (3^12.5).  Hiis  card  is  required  for  on-  and  off-body 
streamlines  as  w^l'll  as  isobars. 

ON- BODY  STREAMLINES  require  no  further  information. 

OFF-BODY  STREAMLINES  require  the  following  card  in  format  (6F12  5): 

DESCRIPTION 

X-coordinate  of  the  dividing  plane  origin. 
Y-coordinate  of  the  dividing  plane  origin. 
Z-coordinate  of  the  dividing  plane  origin. 
X- component  of  the  normal  vector  of  the 
dividing  plane. 

Y-component  of  the  normal  vector  of  the 
dividing  plane. 

Z- component  of  the  normal  vector  of  the 
dividing  plane. 


IAMFTER 

COLUMN 

XI 

1-12 

Y1 

13-24 

11 

25-36 

XN 

37-48 

YN 

49-60 

ZN 

61-72 
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Set  (XN,YN,ZN)  ■  (0.0, 0.0, 0.0)  or  leave  blank  if  all  parts  of  all 
streamlines  are  to  be  plotted. 

ISOBARS  require  the  following  cards.  The  first  is  in  format  (315): 


PARAMETER 

COLUMN 

DESCRIPTION 

NSP 

4-5 

Number  of  body  sections  on  which 
are  to  be  plotted.  (NSP  5  25) 

isobars 

NCPG 

9-10 

Number  of  isobars  to  be  plotted. 
(NCPG  5  25) 

LIM 

15 

If  LIM  f  0,  isobars  will  be  connected 
across  section  boundaries.  Set  LIM=0 
or  leave  blank  to  suppress  this  option. 

The  next  card(s)  is  in  format  (1615)  and  contains  the  section 
numbers  on  which  isobars  are  to  be  plotted. 

The  last  card(s)  is  in  format  (8F10.6)  and  contains  the  values  of 
the  desired  isobars. 


INPUT  CHECKING  AND  PROGRAM  PRINTOUT 


The  printed  output  from  XYZPLOT  lists  input  errors  and  provides  a 
record  of  the  numerical  values  used  to  generate  each  view.  Most  errot 
messages  are  self-explanatory.  Two  of  them  will  be  outlined  here 

1.  THE  LINE  OF  SIGHT  IS  NOT  PERPENDICULAR  TO  THE  PLOTTING  PLANE. 

This  message  is  printed  whenever  the  angle  between  the  line-of-sight  and 
the  xs-axis  differs  from  90  degrees  by  more  than  1  degree.  This  difference 
is  printed  as  "(SKEW-XX  DEGREES)".  A  small  value  of  "SKEW"  will  not  cause 
appreciable  distortion. 

2.  OBSERVATION  POINT  LIES  TOO  CLOSE  TO  BODY.  This  message  will 
occur  during  perspective  projection  when  XYZPLOT  attempts  to  project  a  point 
which  is  very  near  the  observation  point  or  when  the  line-of-sight  :ntet 
sects  the  plotting  plane  at  infinity. 


A  detailed  printout  is  provided  when  plotting  isobars  so  the 
user  may  pinpoint  the  location  of  a  desired  isobar.  Raster  coordinates 
(the  coordinates  used  on  the  actual  picture)  of  the  points  on  each  isobar 
are  listed,  as  well  as  the  body  section  on  which  they  lie.  The  raster 
coordinate  system  is  shown  in  Figure  12. 

Top  of  4060  Output 


Figure  12.  Raster  Coordinates 


COMPUTER  REQUIREMENTS 

XYZPLOT  is  coded  in  CDC  FORTRAN  EXTENDED.  Some  minor  reprogramming 
would  be  required  to  make  it  strictly  ANSI  FORTRAN.  The  program  compiles 
in  approximately  30  CP  seconds  on  the  CDC  6400  processor  and  requires  a 
maximum  of  104K  octal  words  of  storage  to  execute.  If  isobars  are  not 
going  to  be  plotted  in  a  given  run,  0VERLAY(2,S)  (program  ISOBAR)  and  its 
two  subroutines,  CONNECT  and  LOCATE  may  be  omitted,  reducing  the  storage 
requirement  to  75K  octal  words. 
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XYZPLOT  utilizes  a  set  of  plotting  routines  (SCORS  package)  which 
resides  in  both  the  DTNSRDC  CDC-6600  and  CDC-6700  system  libraries.  Three 
scratch  disk  files  are  required  in  addition  to  a  plot  file,  the  three 
possible  data  files,  and  the  regular  input/output  files. 

The  output  for  the  SC-4060  microfilm  plotter  is  written  to  the 
DTNSRDC  batch  plot  queue,  which  is  duiped  to  tape  nightly.  To  obtain  plot 
output,  the  user  must  submit  an  ADP  Off-Line  Request  at  the  ADP  Control 
Center  (the  request  forms  are  available  at  the  ADP  Control  Center  window) . 
Figure  13  shows  an  ADP  Off-Line  Request  with  the  required  information 
shown  by  question  marks.  The  user  will  generally  check  both  the  film  and 
hardcopy  boxes  (in  order  to  obtain  both  film  and  paper  output) ,  and  the 
Tape  ID  must  be  FILMPL-6700  or  FILMPL-6600,  depending  on  which  machine 
the  job  was  run  on.  The  number  of  frames  is  supplied  on  the  XYZPLOT 
printout. 


ADP  OFF. LINE  REQUEST  NDW-HSRDC.  104*2/11 


KEYPUNCH 

□  »  □ 

I  j  vENirv 
!  i  convent  to 

LD  »  □  *• 

□  l-'ST 

I  1  ACAAOOUCt 
,  I  INTeAAAfT 
[  I CNDAUNCM 


OCA 

C  ooc  r*Pe 

[~1  FAINT 

Q] punch 
aoutinC _ 


SC  40*0 
l£  Film 

NO  FRAMES . 


SAtCIAL  INSTRUCTIONS 


CHAN  SIZE  . 

no  ru.es _ 

NO  COA'CS 


OAtAATOA  COMMINTt 


Figure  13.  ADP  Off-Line  Request 
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Table  3  shows  the  DTNSRDC  SCOPE  3.4  control  card  setup  for  a  compile, 
load,  and  go  type  of  nn. 


TABLE  3  -  COMPILE,  LOAD  AND  00  CONTROL  CARDS 


JOBNAME  ,00.04000 ,+  T250. 

•User’s  JOB  CARD* 

CHARGE ,CXXX .PPPPPPPPPP , RS ,B . 

•User’s  CHARGE  CARD* 

RFL.45000. 

•CXXX  is  user's  ID,  PPPPPPPPPP  is  user's 
job  number* 

ATTACH, TAPE2 .BODYPOINTS, ID-CXXX. 

•SURFACE  ELEMENT  POINTS  (if  read  from 

A  disk  file!* 

ATTACH ,TAPE4 .OUTLINE , ID-CXXX . 

•OUTLINE  POINTS  (if  read  from  disk  file)* 

ATTACH, TAPE50 ,nATA, ID-CXXX. 
FTN,T,A, PL-20000. 

ATTACH, SC4020. 

RFL,  104000 .+ 

LDSET , LIB-SC4020 , PRESET-NGINF . 

LGO. 

"  7/8/9  END  OF  RECORD 

XYZPLOT  SOURCE  DECK 

"  7/8/9  END  OF  RECORD 

DATA 

"  6/7/8/9  END  OF  FILE 

•STREAMLINE/ISOBAR  Data  (if  required)* 

75000  if  OVERLAY (2, 5)  and  its  subroutines  are  removed. 


If  the  program  is  to  be  used  often,  it  is  advisable  to  execute  it 
from  an  absolute  overlay  file.  The  following  deck  setup  may  be  used  to 
create  this  file: 
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JOBNAME  ,01104000  ,T50 . 

CHARGE , CXXX , PPPPPPPPPP , CC , B . 

FTN,T,A,PL-10000. 

REQUEST, GOFILE, *PF. 

ATTACH, SC4020. 

LDSET , L I B-SC4020 , PRESET-NG INF . 

LQAD.LGO. 

NOGO. 

CATALOG, GOFILE ,XYZPLOTGOFILE , ID-CXXX .AC-PPPPPPPPPP . 

"  7/8/9  END  OF  RECORD 
XYZPLOT  SOURCE  DECK 
"  6/7/8/9  END  OF  FILE 

Execution  may  then  be  accomplished  by  the  following  two  control  cards: 
ATTACH , GOF I LE , XYZ  PLOTGOF I LE , I D-CXXX . 

GOFILE. 

GOFILE  is  the  local  file  name  of  the  XYZPLOT  overlay  file  and  is 
the  only  name  that  may  be  used. 

Execution  times  are  difficult  to  estimate  because  of  the  large  number 
of  factors  involved.  A  conservative  time  estimate  would  be  50  seconds  per 
view,  plus  30  seconds  for  compilation.  A  table  of  execution  times 
accompanies  the  examples  presented  in  the  next  section. 


INPUT/OUTPUT  EXAMPLES 


Applications  of  XYZPLOT  are  demonstrated  in  this  section.  Output, 
examples  are  given  for  each  of  its  various  uses,  and  the  input  which 
generated  the  on-body  streamline  example  is  shown. 

The  first  exanple  (Figure  14)  contrasts  an  orthographic  projection  witl- 
a  perspective  projection  from  the  same  point-of-view.  This  body  was  a  test 
example  in  a  recent  report  by  Woodward,  et  al.^  The  remaining  examples 
are  in  perspective. 

Woodward,  F.A. ,  Dvorak,  F.A. ,  and  Geller,  E.W.,  "A  Computer  Program  for 
Three-Dimensional  Lifting  Bodies  in  Subsonic  Inviscid  Flow,"  Flow  Research, 
Inc.,  Kent,  Washington,  Rept.  USAAMRDL-TR-74-18  (April  1974). 


Figure  15  shows  a  sphere  with  off-body  streamlines.  The  size  ratio 
(SR  parameter  on  the  sixth  plotting  information  card)  was  set  to  0.4  to 
allow  room  on  the  picture  for  the  streamlines  to  pass  around  the  body. 


Figure  15.  Sphere  with  Off -Body  Streamlines 


On-body  streamlines  are  shown  on  a  sample  helicopter  fuselage  in 
Figure  16.  The  input  is  shown  in  Table  4.  Body  and  streamline  data 
were  read  from  disk  files. 

The  isobar  example  is  Figure  17.  It  shows  an  isobar  on  a  helicopte" 
rotor  hub  fairing:  the  upper  picture  shows  the  results  without  the 
connection  option;  the  lower  picture  shows  a  successful  application  of  the 
connection  option  to  connect  the  isobars  across  the  section  boundaries. 
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Figure  16.  On-Body  Streamlines  on  Sample  Helicopter  Fuselage 


TABLE  4  -  INPUT  FOR  ON-BODY  STREAMLINE  EXAMTLE 


CXXX.CH1Q4003.T50. 

CHARGE. CXXX.PPPPPPPPPP.RS.B. 
ATTACH.GOFILE.XYZPLOTGOFILF  ,IO=CXXX. 
ATTACH,TAPE5J,TAPEl7,IO*CXXX. 
ATTACH,TAPE2,BOOYPOINTS,IO=CXXX. 
SETCORE.INOEF.AOOR. 

GOFILE. 

"  7/8/9  END  OF  RECORD 

3N-B3DY  STREAMLINE  ILLUSTRATION 


•USERS  JOB  CARO* 

•USERS  CHARGE  CARO* 

•XYZPLOT  ABSOLUTE  OVERLAY  FILE* 
• T APE  1 7  FROM  XY?PF»* 

•SURFACE  ELEHEHT  POINTS* 


NANE 


2  0 

CODE 

BLANK  CARO 
1  3 


I 

-20. 

-2C. 

-<*C. 

0. 

9  1 

t  2 

-9. 

7/8/9 
6/ 7/8/9 


2  1 
PHONE 


2C  . 


4  1  576  0. 

(NO  TITLE  ON  FRAME) 

7.3 

0. 

0. 


-1 , 


ENO  OF  RECORO 
ENO  OF  FILF 


(IF  MORE  THAN  ONE  VIEN  NAS  DESIRED,  THE  NEXT  SET  OF  PLOTTING 
INFORMATION  CAROS  MOULD  FOLLOW  THE  LAST  ENO  OF  RECORD  CARO.  ETC.) 
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Figure  17.  Isobar  on  Helicopter  Rotor  Hub  Fairing 


Figure  18  shows  the  use  of  XYZPLOT  in  plotting  body  outlines,  pathlines, 
and  vortex  lines.  These  plots  were  useful  in  visualizing  flow  in  a  rotating 
tank.6  The  tank  was  specified  by  outline  point  and  the  program  which 
generated  the  path  lines  and  vortex  lines  wrote  them  on  a  file  using  the 
XYZPLOT  off -body  streamline  format.  They  were  then  plotted  as  if  they 
were  off-body  streamlines.  Figures  1  and  2  (minus  labeling)  were  also 
generated  using  the  outline  facet  of  the  program. 


TABLE  5  -  LIST  OF  RUNNING  TIMES  FOR  SAMPLE  PROBLEMS 


ILLUSTRATION 

NO.  OF  VIEWS 

EXECUTION 

TIME* 

Figures  162 

2 

19 

Figure  13 

2 

23 

Figure  14 

1 

20 

Figure  15 

1 

40 

Figure  16 

2 

46 

Figure  17 

4 

30 

*  Time  is  in  CDC  6400  central  processor  seconds  (includes  loading  time; 
does  not  include  compilation  time) . 
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FRONT  VIEW 


Figure  18.  Rotating  Tank  Pathline  and  Vortex  lines 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1  DTNSRDC  REPORTS,  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM 
INARY,  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY  CASE 
BASIS. 


